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Atomic and molecular phases through attosecond streaking 
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In attosecond streaking, an electron is released by a short xuv pulse into a strong near infrared 
laser field. When the laser coupling between two states in the target is weak relative to the detuning, 
the streaking technique, which allows for a complete determination of the driving field, also gives an 
accurate measurement of the relative phase of the atomic or molecular ionization matrix elements 
from the two states through the interference from the two channels. The interference may change 
the phase of the photoelectron streaking signal within the envelope of the ir field, an effect to be 
accounted for when reconstructing short pulses from the photoelectron signal and in attosecond 
time-resolved measurements. 

PACS numbers: 32.80.-t, 42. 50. Hz, 42. 65. Re 



With the ability to control and shape infrared (ir) 
laser pulses, a number of new posibilities are emerg- 
ing. The ability to create and characterize intense few- 
cycle phase stabilized pulses [1^ has allowed the formation 
of isolated attosecond pulses [1-3], that in combination 
with the original ir pulse has opened the research field 
of attoscience[^]. In attoscience the quest is for mea- 
surements with attosecond temporal resolution [sl, Q and 
simultaneous angstrom spatial resolution, ideally retriev- 
ing both amplitude and phase of the atomic and molecu- 
lar wavefunctions in real time as an electron moves [slQ- 

The retrieval of atomic and molecular phases requires 
interference between different channels and has recently 
been pursued in the context of two-photon ionization 
with a train of attosecond pulses combined with an ir 
laser pulse [8,, 2\- In these experiments, the measured 
atomic phase is the phase of the two-photon matrix ele- 
ments, which involves a weighted sum over the complete 
spectrum of field- free intermediate states. Other works 
measure multiphoton strong-field ionization phases using 
the process of high-order harmonic generation [10]. In 
this work, we show how to measure the phase of the fun- 
damental one-photon ionization matrix element from an 
initial state to a continuum final state. These atomic or 
molecular phases may serve as valuable tests for detailed 
calculations of the ionization process, and may give direct 
experimental access to the variation of scattering phase 
shifts with energy - a quantity determining the Wigner 
time delay ^^. 

We consider a situation, where the initial state |1) is 
coupled weakly by the ir field to a detuned excited state 
1 2) such that only part of the initial population is trans- 
ferred to 1 2) and only while the driving ir pulse is on, adi- 
abatic following the ir pulse envelope. The relative phase 
between the amplitudes of |1) and |2) will be phase-locked 
to the ir field. We then use an attosecond pulse with a 
duration much shorter than the optical period of the ir 
field to ionize from both the levels, and the two contri- 
bution will overlap in energy and add coherently in the 
overlapping region (Fig. 1). We will show that measur- 
ing the relative phase of the two contributions and at the 



same time the ir field through the attosecond streaking 
process [1] allows a reconstruction of the relative atomic 
or molecular phases from the two ionization channels. 




FIG. 1: (Color online) Schematic view of the ionization pro- 
cess. The two levels |1) and |2) are coupled by the ir laser 
field, and the spectral width of the attosecond pulse is suf- 
ficiently large that photoelectrons from both levels overlap 
energetically. The interference between the two ionization 
channels in the overlapping region carries a signature of the 
phase between the ionization matrix elements. 

The two-level system is described by the wavefunction 
[atomic units (a.u.) with h = e = a^ = rrie = 1 are used 
throughout unless indicated otherwise] 



M/(t) = ci(t)e-^^^*|l)+C2(t)e 



-lE^t 



|2) 



(1) 



with the initial conditions ci(— 00) = 1, C2(— 00) = 0. 
The total Hamiltonian of the system is H = Hq -f VL{t), 
where the two levels are Eigenstates of the field-free 
Hamiltonian Hq and Vl is the interaction of the system 
with the ir laser field. The coefficients Ci{t) may be found 



by solving the set of coupled differential equations 



ci(i) = -i{l\VUt)\2)i 



i{Ei-E2)t 



C2{t) 



(2) 



C2{t) = -i(2|l^i(t)|l)e^(^^-^^)*ci(t). 



Analytical approximations may be applied to these equa- 
tions, giving interesting insights that will be detailed else- 
where. In this work, we parametrize the interaction as 
(1|Vl(^)|2) = pFL{t)^ where p is the coupling strength 
and FL{t) is the electric field strength of the ir field. We 
then solve the equations numerically. 

We are interested in a regime with large detuning, 
S = E2 — El — ujl > pFl^ where ujl is the central laser 
frequency. In this regime, there are no Rabi oscillations 
and the population only follows the field adiabatically. 
After the pulse, virtually all of the population is left 
in the ground state again. When the coupling is much 
smaller than the detuning, the solution to (J2]) may be 
approximated with ci{t) ~ 1 and C2{t) ~ — P^°(^)e'^y25, 
where F^(t) is the envelope of the electric field, such that 
FL{t) = Fl{t) cos{ujLt + (/>). In this regime, the phase- 
evolution of the two terms in ([1]) is determined by the 
laser frequency and the amplitude in the upper level fol- 
lows the envelope Fl{t). 

We now add a second, weaker xuv attosecond pulse 
which ionizes from the state ([!]). The duration of the 
xuv pulse Tx <C cj^^ is sufficiently short that all of the 
electrons will be released at a well-defined phase of the 
ir field. We vary the time, to, of the xuv pulse relative 
to the ir pulse and calculate the electron distribution in 
first order perturbation theory through dP/dkf = iTj^p, 
with 



'-fi 



{to) = -ildt{kf{t)\Vx{t-to)\^{t)), 



(3) 



Here, we approximate the time-dependence of the fi- 
nal state with the Volkov-phase, such that \k{t)) = 
\k)e~i J ^^ ik-\-A{t )) ^ wJiere the vector potential is given 
by FL{t) = —d^/dt and k is the asymptotic momentum of 
the scattering state. In the strong-field approximation \k) 
is approximated by a plane wave, neglecting the Coulomb 
potential in the final state, but this is not necessary for 
the present theory. We factorize the matrix elements as 
{k\Vx{t - to)\n) = Mn{k)Fx{t-to), where n = 1,2 from 
([1]) and Mn(k) is the ionization matrix element for going 
from the state n to the final state |^), which we wish to 
determine. With this factorization and using ([1]) we have 

{k)\Vx(t - to)\^(t)) = En=l,2 Cn{t)Mn{k)Fx{t - to). 

The attosecond xuv pulse has a Gaussian envelope 



<i(i±ii) 



and is given by Fx{t) = Fx^oe ^Ici+e^) g-^^^^^^ where Ux 
the central frequency and ^ is the dimensionless chirp 
of the pulse, which can be measured with the streaking 
technique [12]. The ir pulse has a sine-squared envelope, 

Fhit) = sin^ (7^) cos(cjl^ + 0), which allows the inte- 
grals in the Volkov phase to be calculated analytically. 



As an example, we calculate the ionization spectrum 
in an atom with p = ao and energy difference of 1.85 eV 
corresponding to the energy of the first excited state in 
lithium. The ionization potential of the ground state is 
5.39 eV. For the (near) ir field we choose a 800 nm, 50 fs 
pulse with peak intensity 1 x 10^^ W/cm^. The attosec- 
ond pulse has a full width at half maximum duration of 
290 as and a central frequency corresponding to 91 eV. 
We assume that the electrons are detected parallel to the 
polarization of the ir field. In the first example, the two 
atomic ionization matrix elements are chosen to be 7r/2 
out of phase, such that M2{k) = iMi{k). If the final state 
\k) is assumed to be a plane wave, this would be the rela- 
tion between two states of opposite inversion symmetry 
along the line of detection. Since the two contributions 
are from different initial states, they do not have the 
same symmetry in the ionization continuum. Any phase 
difference would manifest itself as a deviation from the 
^/2 phase difference between the matrix elements. The 
results are presented in Fig. [2l 
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FIG. 2: (Color online) (a) Photoelectron spectrum as a func- 
tion of time delay between the ir and the attosecond pulse, 
from the laser-driven two- level system discussed in the text. 
(b) The center of energy (full, red) and the energy integrated 
ionization yield (dashed, blue). The two curves are oscillating 
in phase as a signature of the ^/2 phase difference between the 
two matrix elements Mi and M2. 

The population transfer with the given set of parame- 
ters is always less than 7%, but still the ionization yield 
oscillates by around 40%. In the limit of weak coupling 
the amplitude of the oscillation is given by p^°(^)/5. In 
Fig. m the population in the excited state is shown along 
with the weak limit estimate. The population left in the 
excited state after the pulse is of the order of 10~^, such 
that the target is left in the ground state. When the laser 



coupling between the two states is smaller than the de- 
tuning, such that there are no Rabi oscillations, the phase 
evolution of the excited level will be locked to the laser 
pulse. This is seen in Fig. [2] (a) by the fact that the max- 
ima in the emission probability coincide with the maxima 
in the electron energy, which follows the vector potential 
of the ir field through the classical electron acceleration 
A^ = - /^^ FL{t)dt = -A{to) and AE{to) - -k ■ A{to). 
Since the energy-integrated ionization probability is de- 
termined by the relative phase of the two contributions 
originating from the ground state and the excited state, 
the ionization probability will be locked to the laser field. 
The advantage of the attosecond streaking technique is, 
that it allows the precise determination of the phase of 
the ir fieldpj, ll3|. Hence, as we can simultanously mea- 
sure the phase of the ir field and the relative phase of the 
electrons ionized from the two levels, we can determine 
the relative phase of the ionization matrix elements Mi 
and M2. 
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FIG. 3: (Color online) The numerically calculated population 
in the excited state, |c2(t)|^ (full, blue) and the approximate 
result from C2(t) ^ —pFL(t)e^ ^ j^s (dashed, red). 

For the same range of parameters as above, but with 
the two ionization matrix elements in phase, M\ = M2, 
the photoelectron spectrum is plotted in Fig. HI As can 
immediately be seen from the photoelectron spectrum, 
the maxima of the ionization probability now coincide 
with half the zero-crossings of the vector potential, show- 
ing that the relative phase between the vector poten- 
tial and the ionization yield is a strong measure of the 
relative phase between the ionization matrix elements. 
This is more clear when looking at the center of energy 
along with the energy integrated ionization probability 
in Figs. [2] and H] (b). In Fig. [21 the two curves oscillate 
in phase, while in Fig. |4]the two curves oscillate ^/2 out 
of phase; the difference in the phase of the oscillations 
maps directly to the difference in the relative phase of 
the ionization matrix elements. 

All of this applies equally well to atoms and molecules, 
as long as there is a coupling between the ground state 
and an excited state. To accurately retrieve the atomic 
or molecular phases, it is necessary to determine the pre- 
cise phase of the vector potential. For polar molecules, 
the center of energy of the photoelectrons does not ex- 
actly follow the vector potential, but may be shifted due 
to the linear Stark shift of the initial state [1^. In this 
case, in order to accurately retrieve the phases of the ion- 
ization matrix elements, either a nonpolar reference state 
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FIG. 4: (Color online) (a) Photoelectron spectrum from the 
laser driven two- level system discussed in the text, similar 
to Fig [21 but with the two matrix elements in phase, (b) 
The center of energy (full, red) and the energy integrated 
ionization yield (dashed, blue). The two curves are oscillating 
f out of phase as a signature of the two matrix elements Mi 
and M2 being in phase. 



or knowledge of the dipole moments is needed in order 
to correct for these effects. 

Effects similar to the polarization effects [l^ may arise 
from the initial state ([1]). If the ionization probabil- 
ity and the vector potential are not exactly in phase, 
that is, when the ionization matrix elements are not ^/2 
out of phase, then the difference in energy of the two 
contributions leads to a shift in the center of energy of 
the electrons. Since the dominant contribution is from 
the ground state, a constructive interference will lead to 
an increase of the average electron energy while a de- 
structive interference will lead to a decrease of the elec- 
tron energy. Now, as the interference is phase-locked to 
the driving laser field, this means that the average elec- 
tron energy will be changed from the classical estimate 
AE{to) = —k-A{to) and the electron energy is not a de- 
finate measure of the phase of the vector potential. This 
may be seen in Fig. [5l where the center of energy is plot- 
ted along with the classical estimate for the same system 
as in Fig. [H but for a laser intensity of 1 x lO^^W/cm^ 
and with a laser wavelength of 750 nm, closer to the res- 
onance. 

In this example a shift of 0.13 rad or 50 as is found. 
We find that in the limit of large detunings, the phase 
shift of the electrons is not sensitive to the intensity of the 
laser. In a simple model for weak couplings to be detailed 
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FIG. 5: (Color online) (a) Photoelectron spectrum from the 
laser driven two level system, parameters are presented in 
the text, (b) Part of the center of energy (full, red) is plot- 
ted along with the classical estimate of the energy variation 
(dashed, black). The two curves are displaced by 50 as or 
0.13 rad 



elsewhere the phase shift in radians may be estimated as 



A = 



28k ■ CL '' 



(4) 



where e^ is the laser polarization vector. For the range of 
parameters used above, this estimates the phase shift to 
be :^ 0.11 rad corresponding to a time shift of A/cj^ = 42 



as. 

That fact that the center of energy of the photoelec- 
trons is not exactly in phase with the vector potential 
may be important when the attosecond streaking tech- 
nique is used both to retrieve atomic or molecular phases 
and when it is used for time-resolved measurements. It 
may give rise to an apparent shift between valence elec- 
trons where near resonances are present and core elec- 
trons with no resonances. The time delay of 50 as in the 
example above is comparable to the temporal resolution 
one can achieve [sl, Q and must therefore be taken into 
consideration. These shifts may be circumvented either 
by calculating the phase shift in the photoelectrons and 
correcting for it or by using photoelectrons from a differ- 
ent level to determine the phase of the vector potential. 
This level should have no near resonances and be nonpo- 
lar to most easily extract the vector potential. 

In scattering theory, a k gradient of the scattering 
phase is known to give rise to a Wigner t ime- delay [ll|, 
which may manifest itself as a time delay in attosecond 
streaking experiments [5]. The Wigner time delay is a 
purely phase-dependent effect and detailed knowledge of 
the phase of the ionization matrix element may lead to a 
better understanding of the importance of this effect. 

In conclusion, we have shown that when an ir pulse 
with a large detuning excites a two-level system, the 
phase evolution of the system will be locked to the pulse 
and the ionization probability will oscillate with the driv- 
ing frequency. The phase of this oscillation is determined 
by the phase of the driving laser field and the phase of the 
atomic or molecular ionization matrix elements. With 
the attosecond streaking technique, it is posible to ex- 
tract the phase of the driving field and hence use these 
oscillations to measure the atomic or molecular phases 
which have so far resisted experimental investigation. 
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